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ABSTRACT 

Contrary to the common view voids have very complex internal structure and 
dynamics. Here we show how the hierarchy of structures in the density field inside 
voids is reflected by a similar hierarchy of structures in the velocity field. Voids defined 
by dense filaments and clusters can de described as simple expanding domains with 
coherent fiows everywhere except at their boundaries. At scales smaller that the void 
radius the velocity field breaks into expanding sub-domains corresponding to sub- 
voids. These sub-domains break into even smaller sub-sub domains at smaller scales 
resulting in a nesting hierarchy of locally expanding domains. 

The ratio between the magnitude of the velocity field responsible for the expansion 
of the void and the velocity field defining the sub voids is approximately one order of 
magnitude. The small-scale components of the velocity field play a minor role in the 
shaping of the voids but they define the local dynamics directly affecting the processes 
of galaxy formation and evolution. 

The super-Hubble expansion inside voids makes them cosmic magnifiers by 
stretching their internal primordial density fiuctuations allowing us to probe the small 
scales in the primordial density field. Voids also act like time machines by "freezing" 
the development of the medium-scale density fiuctuations responsible for the forma- 
tion of the tenuous web of structures seen connecting proto galaxies in computer 
simulations. As a result of this freezing haloes in voids can remain "connected" to 
this tenuous web until the present time. This may have an important effect in the 
formation and evolution of galaxies in voids by providing an efficient gas accretion 
mechanism via coherent low- velocity streams that can keep a steady inffow of matter 
for extended periods of time. 

Key words: Cosmology: large-scale Structure of Universe, methods: data analysis, 
N-body simulations 



1 INTRODUCTION 

Voids are the most salient feature of the Cosmic Web. 
Early galaxy maps unveiled a far from homogeneous Uni- 
verse dominated by vast semi-spherical empty regions sur- 
rounded by a complex network where galaxies aggre- 
gate into de nse massive clusters joined by long filamen- 
tary bridges. (JJoeveer fc Einastdll9 78l:^Tarenghi et al.' 1978; 
Tifft fc Gregorvl Il978l: lEinasto et a l. 1980; Tull v fc Fisher 
19871 : iGeller k Huchralll989l l. This "Cosmic Web" can be 
described as an interconnected network with four basic com- 
ponents: spherical clusters connected by elongated filaments 
defining two-dimensional sheets forming a "net" that sur- 
rounds the voids that account for most of the volume in 
the Universe. This rich network was already imprinted in 
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the primordial density field as tiny Gaussian fluctuations 
according to inflationary models. The development of the 
Cosmic Web from the primordial nearly homogeneous fleld 
is a natural consequence of the gravitational a nisotropic co l- 
lapse as described by the seminar work of Zeld ovichI (|l970h . 
A more explicit description of the relation between the dif- 
ferent elements of the Large S cale Structu re (LSS) came 
with the "bubble" cosmology o f llckd(| 1984 1 and the "Cos- 
mic Web" theory of Bond et al. [1996). These works among 
many others (see van de Weygaert (2002) for an excellent 
review) give insight on particular aspects of the anisotropic 
gravitational collapse that ultimately produces the observed 
cosmic structures but can not fully describe the develop- 
ment of structure into the full non-linear regime. In order to 
understand the non-linear matter evolution one usually re- 
lies on numerical experiments. N-body computer simulations 
have been used to reproduce and study the distribution of 
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matter on large scal es with ^reat success starting from the 
pioneering work of (IDoroshkevich et alj Il980l : iMelott et al.l 
ll983l :lKlvpin & Shandarin 1983) to the current state-of-the- 
art such as the millennium simulation with billions of mass 
particles 



1.1 Studying the underdense Cosmic Web 

Voids are underdense structures that to a first approxima- 
tion have simple structure and dynamics and as such their 
properties are in general considered well understood. Voids 
are expected to have nearly spherical shapes as result of 
their expansion which erases any original sphericity (|Ickd 
1 1984 1). Their internal dynamics can be well described as 
simple expanding domains where 5 oc — Vv indicating that 
the inner regions of the voids expand faster than the outer 
shell s. This gives c haracteristic density and velocity pro- 
files ([Bardeen et al.ll986ll that have been observed i n com - 
puter simulati ons Ivan de Wevgaert fc van Kampen (Il993h: 
Schaap I (|2007h and o bservations (INasonova &: KarachentsevI 
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201 ih . The work of ISheth k van de Wevgaerj (|2Q04l ) rep- 



resented an important step in our understanding of the 
dynamics of voids by dissecting the different scenarios of 
void evolution based in the excursion set analysis (see also 
IParaniape et al.1 (|2012l ll. 

Voids have been the subject to many studies moti- 
vated by their simple structure an d dynamics, the i r low 
density and ap p arently evolution (IZeldovich et al.[ 
Suto et al.1 Il984l: lEinasto et al.1 Il989l: iKirshner et al 
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1993; van de Wevgaert & van Kampen 1993; Szo moru et al 
_1996; Blumenthal et al. 1992 ). The interest in voids led to 
the development of many void finders (see the void com- 
parison project (IColber^ et al.l 20081) ) am ong which the 
watershed void finder (jPlaten et al.l l2007l ) and the Zobov 
(JNev rinck 2008) methods stand out because of their connec- 
tion of local minima with the topopoly of the density field. 



1.2 Galaxies in voids 

Even though voids are by definition empty in some cases 
one can find galaxies at their interior. In fact what is sur- 
prising is that we do not see a large population of low-mass 
galaxies pop ulating voids and tr acing their inner structure 
(JKlvpin et al . 1999; Moor e et all [1999) and that the galax- 
ies we see are b asically representative of the general popu- 
lation (|Peeblesl 12001. ) . Galaxies in voids have been studied 
in de t ail in the Bootes vo id (Pei mbert Sz Torres- Peimbertl 
I1992I : IWeistrop et~al] [1992. ll995l: ISz omoru et al.' Il996l) 
and t he local void ([Karachentsev et al . 2004; Won g et al.l 



I2OO6I : JKreckel et al.ll2Qlll) and the Sloan Digit al Sky Sur- 



vey (|Patiri et al.1 I2OO6I : JCeccarelh et al.1 120081 ). Recently 
new surveys targeted t o voids have been carried out 
with primissing results ( Kr eckel et al.l I2OIII ). The rela- 
tively simplicity of void shapes and dynamics has been ex- 
ploited to propose their use as cosmq l ogical probes with 
tantalizing results ([Padilla et al.1 120051: iLee fc Park! [JOOgI; 
Park k Leell2007l: iBiswas et al.ll2oTQl : IChongchitnan fc Silk! 
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Table 1. Dimensionality of a Cosmic Web element from the neg- 
ative eigenvalues of the deformation tensor. Positive eigenvalues 
indicate collapse while negative eigenvalues expansion. 



1.3 Scope of this w^ork 

This is the first of a series of papers exploring the hierarchi- 
cal properties of the Cosmic Web. In the following sections 
we discuss the emergence and hierarchical character of the 
Cosmic Web as a whole. However, in the present work we 
focus our analysis on voids. Walls, filaments and clusters will 
be treated in following papers. 



THE COLLAPSING AND EXPANDING 
COSMIC WEB 



The seminal work of IZeldovichI (|l970l ) described the evolu- 
tion of a primordial could of matter as a succession of dy- 
namical and geometrical stages where the dimensionality is 
reduced by the gravitational collapse describe by: 



p(x,t) 



1 



(1) 



p (1 - a(t)Ai) (1 - a(t)\2) (1 - a(t)\z) 

where A3 > A2 > Ai. Gravitational collapse occurs when at 
least one of the eigenvalues is positive. When all three eigen- 
values are positive the collapse is complete, first forming fiat 
pancakes then elongated filaments and finally compact dense 
clusters, following the decreasing order of the eigenvectors 
(see Table [1]). 

2.1 The Expanding Cosmic Web:Inverse 
Zel'dovich dynamics 

The Zel'dovich equation can be also seen from the point of 
view of the dimensionality of the different stages in the grav- 
itational collapse given by the negative eigenvalues. Voids, 
walls and filaments correspond to primordial clouds with 
three, two and one negative eigenvalues respectively. This 
strongly suggest that the dynamical character of each mor- 
phology is given by its expanding eigenvectors. In this in- 
verse Zel'dovich picture voids expand in all directions, walls 
expand restricted t o their pla ne and filaments do similarly 
along their length (|lckdll984l ). This local expansion is pos- 
sible even for structures denser than the global mean such 
as walls and filaments because they are locally underdense 
with respect to their surroundings filaments in the case of 
walls and their surrounding dense groups and clusters in the 
case of filaments. We then generalize the idea of a "valley" 
as a locally underdense D^ region with voids, walls and fil- 
aments corresponding to dimension A/" = 3, 2, 1 respectively. 
By accounting for both the positive and negative eigenval- 
ues of the Zel'dovich collapse we have a complete picture 
of a collapsing cloud (positive eigenvalues) and its internal 
dynamics (negative eigenvalues). 
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Figure 1. Hierarchy of structure in the Cosmic Web. The left panel shows the density field on a thin slice across a 64 Mpch~^ simulation 
with ACDM concordance cosmology. The right panels show the highlighted square inside the slice. The density field in the zoomed region 
was computed from a high-resolution resimulation with Npart = 1024"^ centered in the void region. 



2.2 Hierarchical emergence 

What happens in the Zel'dovich collapse when we consider 
a cloud with substructure? This is a more realistic scenario 
than the idealized Zel'dovich collapse since the primordial 
density field is a superposition of waves. The hierarchical 
character of the Cosmic Web is imprinted in the primor- 
dial density field characterized by its power spectrum. At 
the scales relevant for the development of the Cosmic Web 
seen in t he galaxy distributio n (up to a few tens of Mega- 
parsecs) (|Einasto et al.ll2011al l the power spectrum dictates 
that small fluctuations grow first superimposed on large fluc- 
tuations naturally producing a hierarchical scenario. In the 
linear regime all Fourier modes are independent from each 
other and grow independently (this property will be ex- 
ploited in section [3T]) . 

In a hierarchical scenario an over dense cloud will col- 
lapse while internally smaller fluctuations will grow su- 
perimposed to the main collapse. This is the well known 
cloud-in-cloud scen ario that gives rise t o haloes and and 
their su bstructure (Bardeen et al. 1986; B ond et al.l Il99ll : 
iLacev fc Cole 1993; Sheth et al. 2001). On the other side 
of the Zel'dovich equation we have several scenarios for 
the development of underdense regions as described by 
ISheth &: van de WevgaertI (|2004l l for the case of D^ valleys. 
Here we have a hierarchy of voids where a given underdense 
region can contain locally under dense sub-regions or even 
be contained by a larger overdense region. In the first case 
the void expands while internal sub- voids also expand super- 
posed to their parent voids's expansion. In the former case a 
void can be embedded inside larger overdense systems that 
eventually lead to the collapse of the void itself. 



plex picture that is not seen in th e galaxy distribution 
(iMathis fc Whitell2002l : iGottlober eTal. 2003 : Colberg et al.1 
|2005[ ) . If galaxies populate the dark matter haloes found in 
computer simulations then one would expect to see the usual 
network of galaxies traced by large haloes and underneath 
it a much richer sub- web of tenuous fila ments containing a 
large population of low-mass haloes (Dubi nski et al .11 19931 ). 
The large-scale fluctuations that give origin to cosmo- 
logical voids contain also smaller fluctuations which grow 
inside what behaves like a locally low-density universe. The 
superimposed embedded fluctuations grow faster than their 
parent void although damped by the parent void expansion. 
As a result of this a given D^ valley will contain and be de- 
fined by smaller D^ valleys. This process occurs at all levels 
in the hierarchy of voids as sub- voids themselves also contain 
smaller fluctuations (see Figure [1]). 



2.3.1 Hierarchy of walls and filaments 

The substructure scenarios described above can be difficult 
to visualize in the case of walls and filaments if one only fo- 
cuses on their collapse. However, if one focuses on the nega- 
tive eigenvalues then we can see that an expanding D^ val- 
ley can contain embedded D^ and D^~^ sub- valleys which 
expand locally with respect to their parent expansion. This 
restriction in the dimensionality of substructure comes di- 
rectly from equation [T] and means that a filament can have 
sub-filaments and nodes but no sub- walls or sub- voids. A 
wall can contain sub-walls and sub-filaments but no sub- 
voids and a void can contain all kinds of substructure (see 
Table [2|). 



2.3 The hierarchy of voids 

The distribution of matter inside voids studied with 
high-resolution computer simulations shows a very com- 



2.4 Hierarchical velocities in voids 

The velocity field in the Universe is dominated by large 
streaming or bulk motion. Our own Mily Way together 
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Table 2. Allowed substructures according to the dimensionality 
of the parent structure (D^). Voids can contain all sub-elements, 
walls can contain sub-walls and sub-filaments and filaments can 
contain only sub-filaments. 



with the local group is moving with a speed o f ^ 600 

Dressier! 



km/s towards the great attractor (Shava 198 - ._ 

Il987l : iLvnden-Bell et alJ FlOSSl : lAaronson et alJ Il989l l and 
at 260 km/s in th e opposite direction from the local void 
(|Tu11v et al.ll2008l V At medium scales he dynamics around 
the Milky Way are affected by the local expansion of our 
parent wall, producing the so-called "cold Hubble flow" 
(Arago n-Calvo et al.ll201lh . The power spectrum of the ve- 
locity field Pv{k) is given in terms of the power spectrum of 
the density field Pm{k) as: 



Py{k) (Xk-^Pm{k). 



(2) 



The k~^ boost the power at large scales while the small 
scales have very little power in comparison. From this we 
see that the velocity field is a very good indicator of the 
large scale matter distribution. The velocity fluctuations at 
a given k will be superimposed on much larger fluctuations 
in the form of bulk- flows which originate from larger density 
fluctuations. In fact the velocity field correlates poorly with 
the observed features in the cosmic web (see top panels in 
Figure (8]). The k~^ factor effectively produces a mismatch 
between the dominant scales in the density field and the 
velocity field. 

Consider a filter that applied to the velocity field will 
attenuate the power at large scales such as W oc /c^ in which 
case Pv{k) oc Pm{k). Now the fluctuations in the velocity 
field have the same scaling as the fluctuations in the density 
field and this will be reflected in a closer match between the 
two. Another option is a Gaussian high-pass filter. This filter 
attenuates the amplitude of Pv{k) at large scales bringing 
the shape of the power spectrum closer to Pm{k) after the 
cut-off scale. A high-pass filter provides a way to select scales 
in the velocity field associated to particular elements in the 
Cosmic Web. The high-pass filtered velocity field contains 
only local fluctuations at the cut-off scale. This is analogous 
to measuring peculiar velocities after removing the averaged 
local velocity around a window of a given size. Since the 
high-pass filter matches the velocity and mass power spec- 
trum we will be able to identify small scale features in the 
density field reflected in the velocity field. We should then 
expect to see a hierarchy of expansion domains at the inte- 
rior of voids reflecting the hierarchy of voids observed in the 
matter distribution. 

Equation [2] is valid in the linear regime. At the present 
time small-range gravitational interactions make the pecu- 
liar velocities in the non-linear, high-density environments 
very high and turbulent as power is transfered from large 



scales. On the other hand, the low-density nature of cosmo- 
logical voids results in a velocity field dominated by large 
laminar motions and at smaller scales low- velocity coherent 
streams. This makes voids ideal laboratories to study the 
velocity field at all scales. 



2.4' i Residual velocities 

In order to isolate the local component of the velocity field 
from the bulk flows that dominate the overall dynamics of 
the Cosmic Web we define the residual velocity at scale r as: 



sir) 



VG(r) 



(3) 



where v is the original velocity field and VG{r) is the velocity 
field smoothed with a Gaussian filter of radius r. This is a 
high-pass filter that suppresses flows at scales larger than 
r and allows us to probe the velocity field at an specific 
local scale. In the following sections we will use this filter to 
explicitly expose the hierarchy of velocities in voids. 



2.5 Voids as cosmic microscopes and time 
machines 

Voids offer a unique opportunity to probe the small scales 
in the primordial power spectrum and the early develop- 
ment of medium-scale fluctuations. The small-scale density 
fluctuations embedded in voids are stretched as the void 
expands, effectively turning voids into cosmic microscopes. 
Inside these low-density environments the growth of primor- 
dial fluctuations is attenuated and continue in the linear and 
semi-linear regimes for longer time than in denser regions. 
Voids then offer a direct way to study the primordial power 
spectrum at smaller scales than in their complementary high 
density regions where the opposite effect occurs and small 
fluctuations are compressed and even erased by the gravita- 
tional collapse. 

There is a perhaps more interesting property of the sub- 
structure inside voids originating from the break in symme- 
try in the collapse of small fluctuations depending on their 
position with resp ect to larger fluctua t ions. F ollowing along 
the sam e lines as ISuhhonenko et aU (|201lh : lEinasto et alJ 
(|2011al lb[) lets consider the evolution of a given set of small- 
scale fluctuations in different large-scale environments. If the 
small-scale fluctuations are located on top of a large over- 
dense peak they will begin to grow enhanced by their po- 
sition inside the peak and at the same time defining a full 
medium-scale cosmic web connec ting the small-scale peaks 
(see figure 3 in iRomano-Diaz et al.l (|201l|)). The larger vol- 
ume of the peak will eventually collapse erasing any memory 
of the primordial fluctuations and their medium-scale web. 
On the other hand, if the same small-scale fluctuations are 
located on a large underdense valley their growth will be 
attenuated by the expansion of the void but more impor- 
tant, the expansion will stretch the fluctuations and they will 
never collapse as a whole as in the case of the large peak. 
Instead they will evolve into a tenuous web with its own 
network of voids, walls, filaments and low- mass haloes. The 
expansion of the void will further stretch this inner web to 
larger scales (see Figure [1]). The end result is that the band 
of k in the power spectrum responsible of producing the ten- 
uous cosmic web inside voids is effectively shifted towards 
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higher values of k. Voids then offer an unique opportunity 
to study the early evolution of the "frozen" cosmic web at 
medium- scales. The same idea can be applied to walls and 
perhaps even to less extent to filaments, although the fluc- 
tuations there correspond to the power spectrum collapsed 
along one and two directions respectively, making the inter- 
pretation not as direct as in the case of voids. 



Adjacency condition Structure 

2 voids wall 

^ 2 walls filament 

^ 2 filaments node 



Table 3. Criteria used to classify structures in the Cosmic Spine 
method. 



3 DECONSTRUCTING THE UNIVERSE: THE 
COSMIC SPINE 

The Hierarchical Cosmic Spine is the most recent 
implementation of the origi nal method introduced i n 
lAragon-Calvo et al.1 (|2010l l and lAragon-Calvo et aJ] (|201(t l. 
The main improvements presented here are the addition of 
explicit hierarchical relations for voids, walls and filaments 
and a self-consistent criteria for the identification of cosmic 
structures. 

We start by identifying the voids in the density field us- 
ing a floatin g-point imple mentation of the original watershed 
void finder ( Platen et al.l l2007V This provide us with a par- 
tition of space into contiguous regions sharin g the same local 
minim a. In the Cosmic Spine method ( Aragon-Calvo et al.l 
|2010|) we assign the boundaries between the watershed re- 
gions (i.e. the watershed transform) to the network of walls, 
filaments and clusters. 

From the topology-based segmentation provided by the 
Watershed Void Finder we characterize the local geometry of 
the density field as follows: We start by labeling regions iden- 
tified with the watershed transform as cosmological voids. 
Each void has its own label provided by the watershed trans- 
form. We then focus on the boundaries between void regions 
and label the intersection between two voids as walls. Then 
the intersection between walls correspond to filaments and 
the intersection between filaments are the nodes of the net- 
work which in the Universe correspond to galaxies, groups 
and clusters (see Table [3|. This is done in practice by iterat- 
ing through all the voxels labeled as watershed boundaries 
and then labeling as walls those voxels that have two dif- 
ferent voids in their adjacent 3x3x3 neighborhood. We 
run again over the remaining unlabeled boundary pixels and 
identify filaments as voxels have at least two walls in their 
adjacent neighborhood. Finally, we identity haloes as vox- 
els having two or more adjacent filaments. Even though our 
method provide us with haloes as the nodes of the network 
we opt to identify them using other methods like FoF that 
are not limited by the grid resolution. The Cosmic Spine 
method provides with a complete framework for the charac- 
terization of the cosmic web into its primary constituents: 
voids, walls, filaments and clusters. It has no free parameters 
(with the possible exception of grid resolution which is not 
an intrinsic free parameter since its variation does not affect 
the final result as long as the spacial sampling can recover 
the features in the density field) and being fundamentally 
a topological measure it is highly robust agains noise and 
artifacts in the density field. 



3.0.1 Hierarchical SpineWeb 

The main limitation of the original Cosmic Spine method 
is the lack of control in the structures that we identify. A 



phenomenon referred to as oversegmentation. Being a topo- 
logical method any feature in the density field that aflPects 
its connectivity will be detected regardless its contrast. For 
instance, consider the 2D case of a void with very dense 
boundaries but an extremely tenuous bridge of matter across 
its center. The Spine method will identify two voids as long 
as there is a bridge of matter dividing the slope lines ema- 
nating from the local minima. Several approaches have been 
proposed to deal with these substructures inside voids all in- 
volving the use of a free p arameter as some form of threshol d 
or significance measure (Platen et al.l 120071 : ISousbidl201lh . 
We take different approach and first consider the nature of 
the substructure. In N-body computer simulations where we 
have complete control over the sampling of the underlying 
density field (particles in this case) we expect that the fi- 
delity of the features in the density field is limited by the 
accuracy of the method we use to estimate the density and 
that (beyond numerical accuracy) any feature in the den- 
sity field is the direct result of gravity. From the discussion 
in Section 12.21 we see that substructure is a natural con- 
sequence of the hierarchical development of structure and 
therefore we should not remove it by using ah-hoc proce- 
dures. Instead the proper way to deal with such hierarchy 
of structures is by using a hierarchical approach to the iden- 
tification of structures. This is the main motivation behind 
the Hierarchical Spine Web. By doing a hierarchical classi- 
fication we avoid the need of an extra free parameter and 
instead consider the Cosmic Web as a network of nesting 
structures where large structures contain and are defined by 
smaller ones. 



3.1 Constructing hierarchical Spaces 

The crucial point in the hierarchical analysis is the se- 
lection of a proper filter that can expose the hierarchy 
of structures in the Cosmic Web. This is accomplished 
by doing the Cosmi c Spine analysis in hierarchical-space 
(I Aragon-Calvo et al.ll2010. ) which is a generalization of scale- 
spaces that are defined by one of more properties that are 
manifestations of the hierarchical nature of the Cosmic Web. 
Hierarchical spaces are not restricted to linear scale-space 
operators such as the Gaussian kernel (|Florackl 1 19931 ) . In 
the case of N-body simulations we can take advantage of 
having access to the full evolution of the simulation to de- 
sign a filter that targets specific scales in the linear regime at 
early times when the Fourier modes of the density field are 
independent and grow independently (see Section 12. 2|) . This 
linear-regime low-pass filtered density field will evolve into 
a Universe with all the structures above the cut-off scale in 
place and shaped by the anisotropic gravitational collapse. 
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Figure 2. Density field across a thin sfice of the simulation box at 2; = after linear-regime smoothing for three cases: original field (left 
panel), 2 h~^M^pc (central panel) and 4 h~^M^pc (right panel). 



but lacking of the small-scale details. This approach is fun- 
damentally different from the usual a posteriori smoothing 
operation, in that it avoids the nonlinear effects resulting 
from cross talk between Fourier modes. It has the advantage 
of cleanly exposing the hierarchy of structures imprinted in 
the initial density field. The linear-regime filter has a differ- 
ent impact on each element of the Cosmic Web. In the case 
of voids the initial fluctuations will be smoothed at a given 
comoving scale in the linear regime but as the void evolves 
and expands the embedded fluctuations will be stretched, 
effectively increasing the scale of the smoothing length. In 
the case of clusters we have the opposite effect with the mat- 
ter inside the turnaround radius collapsing and erasing the 
small scale fluctuations. Filaments and walls are intermedi- 
ate cases. 

Figure [2] shows the density field at z = evolved from 
the original initial conditions (left panel) and linear-regime 
smoothed versions at 1 /i~^Mpc and 2 /i~^Mpc (center and 
right panels respectively). The linear-regime smoothed ver- 
sions lack of small scale details. Even though there are no 
small scale features the large voids still have well defined 
sharp boundaries in contrast to what we would see in a 
Gaussian smoothing. At the top of the hierarchy we only 
have large scale structures and basically featureless voids. 
The middle level shows some features appearing inside the 
voids while the bottom level shows a rich inner-web at the 
interior of the voids. This progressive addition of features is 
exploited in the following sections to identify and character- 
ize the full hierarchy of structures in the cosmic web. 



3.2 Hierarchical void finding 

We start the hierarchical characterization of the Cosmic Web 
by identifying voids regions as described in the previous sec- 
tion. In the Cosmic Spine hierarchical framework we identify 
voids independently at all levels in hierarchical space, and 
establish the cross-scale relations between voids at different 
levels. A given parent void at the hierarchy level i is defined 
by smaller children voids at the next hierarchical level, i+ 1. 
We assign parent- child relations by identifying overlapping 
volumes between voids in adjacent levels in the hierarchi- 
cal space. A given child usually shares volume with several 



parent voids higher in the hierarchy. We enforce a non-loop 
property in the hierarchical tree by assigning each child void 
exclusively to the one parent void to whom the child con- 
tributes most of its volume. This constraint assures that all 
children voids have only one single parent in the void tree 
hierarchy. 

We perform the void merging across adjacent levels in 
the hierarchy by computing only the flooding procedure on 
the watershed method (see Aragon-Calvo et al. (2010) for 
details). This procedure segments the density field into wa- 
tershed basins but does not explicitly provide the boundaries 
between adjacent watershed regions. The "incomplete wa- 
tershed" focuses on the space partitioning aspect of the wa- 
tershed transform. This makes it straightforward to merge 
voxels between children and parent voids producing com- 
plete hierarchical void tree. The parent void is then recon- 
structed as the union of its children voids as: 



Vp 



u- 



(4) 



where Vp is the set of voxels defining a parent void and Vi 
are the voxels defining each void child contained by Vp . Note 
that the basins of the incomplete watershed also contain the 
watershed transform. After the merging procedure we com- 
pute the full watershed transform (i.e. watershed basins and 
boundaries) using the reconstructed parent void basins given 
by eq.|3]and the density field with the highest resolution, cor- 
responding to the bottom level of the hierarchy. This new 
watershed contains only large structures corresponding to 
the parent voids while conserving the small scale features 
present in the high resolution density field. The particular 
watershed implementation we use is based on a local flood- 
ing restricted to the voxels adjacent to void boundaries. This 
reduces the number of voxels needed to sort in the flooding 
step and greatly improves performance. We can fully ana- 
lyze a set of three 512^ voxel grids in a few minutes using a 
standard workstation. 



3.3 Hierarchical Cosmic Spine 

From the hierarchy of voids described in Section [3. 2 1 we pro- 
ceed to identify walls, and filaments using the criteria out- 
lined in table [3] for each hierarchical level independently. 
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Figure 3. Hierarchical Cosmic Spine superimposed to its corresponding density field for the bottom, middle and top levels (left, central 
and right panels respectively). The black lines delineate the boundaries of the voids found with the Hierarchical Cosmic Spine. The 
background corresponds to the density field across a thin slice of the simulation box at z = after linear-regime smoothing for three 
cases: original field (left panel), 2 h~^Mpc (central panel) and 4 h~^Mpc (right panel). 



Name 


Npart 


Effective scale 


lOOMpc-AO 


5123 


0.25 Mpc/h (Nyquist) 


lOOMpc-Al 


5123 


1 Mpc/h 


100Mpc-A2 


5123 


2 Mpc/h 


100Mpc-A4 


5123 


4 Mpc/h 


lOOMpc-BO 


2563 


0.4 Mpc/h (Nyquist) 


lOOMpc-CO 


1283 


0.4 Mpc/h (Nyquist) 


lOOMpc-DO 


643 


0.8 Mpc/h (Nyquist) 



Table 4. Effective scales of the simulations used in this work. The 
effective scale indicates the cut-off scale given by the linear regime 
Gaussian smoothing or the natural scale given by the particle 
sampling. The suffix in the name indicates initial conditions 
where the particle sampling corresponds to the Nyquist limit. 



The hierarchical relationship between elements of the Cos- 
mic Web in adjacent hierarchical levels is implicitly given 
by the reconstruction of parent voids from their children 
voids so there is no need to explicitly establish parent-child 
relations again for walls and filaments. 

This hierarchical identification of structures results in 
a watertight network of voids surrounded by walls and fil- 
aments. Each void contains an inner network of sub- voids 
each with its own surrounding system of sub- walls and sub- 
filaments. A more detailed description of this procedure will 
be given in the following papers of this series. In the present 
work we focus only on the identification of voids. 



4 SIMULATIONS 

The analysis presented here is based on an N-body simula- 
tion with dark matter corresponding to the standard cosmol- 
ogy with cosmological parameters Qm =0.3, Qa = 0.7, /i = 
0.73, (78 = 0.8. Initial conditions were generated for 512^ 
particles in a periodic box of 100 Mpc/i~^ using the pub- 
licly available cosmics software (JBertschingerl 1 19951 ) with 
power spectrum computed from an analytical fit. For practi- 



cal purposes like fast visualization and preeliminar analysis 
we generated lower resolution versions at 256^, 128"^ and 
64^ particles followin g the averaging procedure described in 
iKlvDJnet aD(|200lh . 



4.1 Linear-regime smoothing 

In order to compute the hierarchical Cosmic Spine we gen- 
erated linear-regime smoothed versions of the initial condi- 
tions as described in Section 13. 1[ We take a diflFerent ap - 
proach to the one presented in Arag on-Calvo et a l.' (2010y 
were we kept the particle sampling at the Nyquist frequency 
by merging adjacent particles and effectively lowering the 
particle count by a factor of 2 on each dimension at each 
smoothing step. In the present analysis we perform the 
linear-regime smoothing on the original initial conditions us- 
ing a Gaussian window while keeping the particle number 
constant. In practice we assign positions to particles on a 
regular grid using the smoothed displacement field as given 
in the Zel'dovich approximation (j Zeldovich 1970,) : 



x = q + *^(q,2;). 



(5) 



Where ^^(q, ^) is the perturbation field smoothed with a 
Gaussian function at scale a. Velocities are assigned accord- 
ing to: 



v = f{z)n{z)^^{q,z). 



(6) 



Where f{z) is the logarithmic derivative of the growth factor 
and Ti is the reduced Hubble parameter. With the above pre- 
scription we generated linear-regime smoothed initial condi- 
tions at 1, 2 and 4 h~^Mpc (see Table IJ]). We then followed 
the evolution of each of the initial condition boxes from 
z = 49 until t he present tini e, z = 0, using the GADGET-2 
N-body code (|SDringelll20"05h (see Figure [2]). 



4.2 Computing Density field 

From the final particle distribution we compute the den- 
sity field for the three simulations described in the previous 
section. This is done usind the Delaunay Tessellation Field 
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Figure 4. Density field across a thin slice of the simulation box 
(top panel) and its corresponding void network (bottom panel). 
Only the boundaries of the voids are shown as semi-transparent 
membranes each with a different randomly assigned color. We 
show the full tridimensional extent of all voids that are inter- 
sected by the density slice on the top panel. The boundaries were 
visualized by volume rendering the pixels in the boundaries of the 
voids with a random color transfer function. 



Interpolator (DTFE) method ([Schaap fc van de Weygaert 
I2QQOI : ISchaap ll2QQ7l : Ivan de Weygaert fc Scha'ani2QQ9 ) .The 
DTFE method is a parameter- free density estimator and in- 
terpolator. It assigns densities to particles acoording to the 
inverse of their adjacent Voronoi cell as: 



where F (Wi) is the volume of the adjacent Voronoi cell and 
m is the mass of the particle. The ability of the DTFE 
method to follow the intrinsic anisotropics in the density 
field is crucial for the correct identification of structures by 
the Cosmic Spine. 



4.3 Computing Velocity fields 

We interpolate the velocity field into a regu- 
lar grid by means of th e Dela u nay tess el lation 
(JBernardeau fc van de WevgaertI Il996l : ISchaap I l2007l : 
iRomano-Diaz k, van de WevgaertI 120071 ) . The velocities at 
each particle in the simulation are assigned to its corre- 
sponding vertex in the Delaunay tessellation of the particle 
distribution. The velocities are then interpolated from the 
D -\-l Delaunay vertices (per tetrahedron) ro, ri, r2...rAr at 
the central position r of every pixel on a regular grid as: 



v{y) = v{yq) + V^;|Dei • (r - ro) 



(8) 



where \/v\^e\ is the field gradient inside the Delaunay tetra- 
hedron. 

In order to avoid aliasing from multistreaming when 
sampling velocities on a regular grid we smooth the velocity 
field per particle at the corresponding voxel size (Nyquist 
limit) by performing a weighted mean as: 



Er 



WiVi 



En 



where the weights Wi are computed from 
tion: 



(9) 



Gaussian func- 



^/2t\ 



(10) 



p(x) = 



(D + l)m 

viyVi) 



(7) 



and n is the distance to particle i. The weighted averaging 
is done separately per velocity component. This smoothed 
velocity field is then used to interpolate on a regular grid 
using the Delaunay tessellation as described above. 



5 THE HIERARCHY OF VOIDS 

We perform the hierarchical Cosmic Spine decomposition 
using the initial conditions with 0.2 /i~^Mpc (Nyquist) 
and with linear-regime smoothing at 2 and 4 /i~^Mpc 
corresponding to simulations lOOMpc-AO, lOOMpc-Al and 
100Mpc-A2 (see table |4]). Ideally one would want to cre- 
ate a continuous hierarchical space. In practice, however 
this is not possible and one must decide the number and 
spacing of hierarchical levels beforehand. The design of 
scale/hierarchical spaces is a non-trivial task and in general 
the scale-space resolution is limited by the available compu- 
tational resources fsee lSato, Y. et al.l (|1998) for an extended 
discussion) . In our case the creation of the hierarchical space 
involves the full evolution of a 512^ particles simulation 
which limits the number of levels that one can produce and 
analyze in a reasonable amount of time. In our experience 
three hierarchical levels give a good compromise between 
scale coverage and the computational resources needed. Note 
that in the case of voids the dynamic range of the attenu- 
ated scales covered by the hierarchical space is larger than 
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Figure 5. Hierarchy of sub- voids inside a parent void. The top panel shows a shce across the density field from which the central region 
with void is highlighted with a write square. The central void is then decomposed into its embedded sub-voids in the middle panels. The 
tenuous gray surface corresponds to the boundaries of the parent void while the colored surfaces to the sub- voids. We show three sets if 
sub-voids with the right middle panel showing all the sub-voids contained in the parent void. Note that the parent void contains and is 
defined by its sub- voids. Bottom panels focus on the largest sub- void (marked with a letter A in the middle right panel) and decomposed 
it on its sub-sub- voids in a similar way as in the middle panels. 



the initial range imposed by the linear-regime smoothing be- 
cause the scale of the density fluctuations is amplified inside 
voids giving a larger effective smoothing length. 

In what follows we denote the hierarchical levels ob- 
tained from the original density field and the linear-regime 
smoothed field at 2 and 4 h~^Mpc as top, middle and bot- 
tom levels respectively. 

Figure |4] shows the void network identified at the top of 



the hierarchy. The top panel shows a slice through the den- 
sity field and the bottom panel shows the network of voids 
that cross the slice. For clarity we show each void with a 
different color. The void network covers all the space in the 
Cosmic Web. The net of clusters, filaments and walls lies 
at the boundaries of the void network. Voids have complex 
shapes, not only departing from the ideal sphericity but also 
their boundaries have a very complex "texture" reflecting 
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the non-linear nature of the dense boundaries. It may be 
surprising that voids are not regu lar spher es as one would 
expect for an expanding region (|Ickelll984l ). However, this 
is valid only for an isolated void. In the context of the Cos- 
mic Web the shape of an expanding void is affected by the 
expansion of its adjacent voids and even more complex pro- 
cesses like the collapse of the void itself (the void-in-cloud 
scenario) in which case the shape of t he void would be highly 
elongated (Sheth & van de We vgaeT^IJOoi l. 

The hierarchical decomposition of a single void is illus- 
trated in Figure [5] The top panel shows the density field 
across a thin slice of the simulation box. A white box high- 
lights an interest region with a void at its center. The cen- 
tral void has well defined dense boundaries and two large 
clusters almost at opposite directions. This particular void 
has significant substructure across the slice. In order to il- 
lustrate the hierarchical decomposition of voids we select 
the central void in Figure [5] (top panel) and show the sub- 
voids it contains by successively adding them to their parent 
void (middle panels). The boundaries of the parent void are 
shown as a thin gray surface and the sub-voids are shown 
as colored surfaces highlighting their boundaries. We then 
take the largest sub-void (marked as A) and decompose it 
in its sub-sub- voids (bottom panels). Here the sub- void acts 
as a parent void in this hierarchical level and is shown as a 
thin gray surface. We add sub-sub- voids as colored surfaces 
just as in the medium panels. There is a large variation in 
the relative size of the sub- void with respect to their parent 
void. The selected void has three large sub- voids (purple, 
green and red) and several smaller sub- voids. The relative 
size of sub-sub- voids seems to be more regular. 



5.1 Void size distribution 

The distribution of void sizes is shown in Figure [6l We es- 
timate the void radius from its total volume as Hvoid = 
(3/(47r)y)^/^. Where V is the volume of ah the voxels in 
the void. The distribution of void sizes is shown for the 
three hierarchical levels computed in this work. The three 
curves have similar shape and are close to a log-normal dis- 
tribution. The peak of the distribution increases as we move 
across the void hierarchy at ^^ 4, 6 and 8 /i~^Mpc for bot- 
tom, middle and top levels respectively. The peak of the top 
level voids is ^ 11 Mpc for h — 0.73. Recently estimates 
of the radius of the local void by Nasonova fc Karachentsev 
(| 20 111 ) give Rlv — 10 Mpc. According to our results the lo- 
cal void is a fairly common void in terms of its size. The 
largest voids we identify in our simulation correspond to the 
top of the hierarchy and have a maximum radius close to 
20 /i~^Mpc. The medium and bottom levels have maximum 
cut-off radius of 12 and 9 /i~^Mpc respectively while we can 
find very small voids in all levels of the hierarchy. The top 
lev el voids in Fi gure [6] are equivalent to the voids identified 
bv lPlaten et al.l ((20071) and in fact our results are in excel- 
lent agreement with their single- level analysis. The top of 
the hierarchy has a very strong correlation with the voids 
delineated by the dense filaments and clusters of galaxies. 
Combini ng; our results with the exc ursion set formalism de- 
rived bv iFurlanetto fc PiranI (|200Q l we expect the voids in 
the three hierarchical levels to be delineated by galaxies with 
a range in luminosity of Mr < — 16 and Mr < — 20 between 
the bottom and top levels respectively. 
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Figure 6. Probability distribution function of void radius for the 
top, middle and bottom hierarchical levels computed in this work 
corresponding to solid, dashed and dotted lines. 
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Figure 7. Density profile of voids as function of enclosing radius. 
The light-gray lines show several individual void profiles while the 
black line corresponds to the mean of all the profiles. The x-axis of 
the individual profiles are normalized with the peak of the curve.. 



5.2 Density profile of voids 

Figure [71 shows a first impression of the internal structure of 
voids. Here the density profile was computed inside increas- 
ing spheres placed at the center of voids. The void centers 
were computed as the mean of all the voxels defining the 
void. The x-axis is scaled such that the density profile for 
each individual void reaches its peak at R=l for easy com- 
parison between different voids. This also allows us to stack 
several void (gray lines) and compute its mean profile (black 
line). The density profile starts around ^ + 1 ^ 0.1 at the 
center of the voids and then increases until it reaches the en- 
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closing network of walls, filaments and clusters at the peak 
of the density profile. After this point the profile slowly con- 
verges to the mean density. Note that some profiles have 
their peak close to the mean density and some even don't 
seem not reach it. These are shallow voids that are embedded 
in very large underdensities of the order of ~ 100 h~^Mpc 
and structures that are still in the process of form ation 
(|Einasto et al.' 1997; She th &: van d e Wevgaert 2004). The 
density profile inside voids is the result of the draining of 
matter from the void into its dense boundaries as they col- 
lapse and increase their density high above the mean. The 
void profiles we com puted are in go od agree ment with theo- 
retica l expectations iHoffman et al.l (|19 83): Hausma n et al.l 
(Il983l): iFillmore &: GoldreichI (Il984l ): fBertschinger. (j 19851 ): 
ISheth fc van de WevgaertI (|2004l l and the density profile re- 
cently computed from the distribution of galaxies in the local 
void by tNasonova fc Karachents ey (.2011, ). 



6 VOID DYNAMICS 

We now focus on the internal dynamics of voids and their 
hierarchical character. Figure [8] shows a slice of the density 
field across the simulation box (top left panel) and the corre- 
sponding velocity field on the remaining panels. The velocity 
field was visualized using the particle advection technique in 
which test particles are randomly placed on the simulation 
box and their positions are evolved following the instanta- 
neous velocity field. The velocity fields shown here corre- 
spond to the two-dimensional velocity field along the x — y 
plane centered on the slice under consideration. The struc- 
tures in this two-dimensional field are easier to visualize than 
in the three-dimensional case where one has structures along 
the line of sight which can lead to confusion. The advection 
field was created by iteratively displacing the test particles 
according to: 



x-^^ =x^ +v(x^)^ 



(11) 



where the index j indicates the step in the advection process, 
x;^ is the position of the particle i and v(x:^):^ the velocity 
field evaluated at the current particle's position. Each par- 
ticle is also assigned a color based on their advection step j 
according to a given lookup table. 

From the slice in Figure [8] we select a region (yellow 
square) with a large void and focus on it. The correspond- 
ing unfiltered velocity field is shown in the right top panel. 
The structures in the velocity field have some correspon- 
dence to the structures see in the density field. There is a 
weak correlation between the centers in the velocity field 
and the center of the large voids although the boundaries of 
the voids are not delineated by the velocity field. The voids 
themselves seem to be embedded in a much larger bulk fiow 
going from the top right corner to the bottom-left corner 
of the box. This is a consequence of the velocity field being 
dominated by fiuct nations at larger scales. The correspond- 
ing amplitude of the velocity field is shown in Figure [9] (left 
panel) . These large scale streams have amplitudes up to 600 
km/s near massive structures. The bulk flow inside the cen- 
tral void of Figure [8] has velocities in excess of 300 km/s. 
There is a weak correlation between the magnitude of the 
velocity field and the underlying distribution of matter. 



6.1 Voids as hierarchically expanding domains 

In order to highlight the contribution to the velocity field 
from density fiuctuations at scales relevant to voids we ap- 
ply a high pass filter to the velocity field (see Section [2.4. II 
for details). We choose filtering scales that fall below the 
peak in the void size distribution at r < 10 h~^ Mpc (see 
Figure (6]) in order to probe the internal dynamics of voids. 
Figure [8] shows the residual velocity field at scales 4 and 1 
h~^ Mpc (right middle and bottom panels respectively). The 
structures in the residual velocity field now show a striking 
correspondence to the structures in the density field. The left 
middle panel shows the central void as a locally expanding 
domain with well defined boundaries closely matching their 
counterpart in the density field. Figure [8] shows voids more 
clearly identified in the residual velocity field than in the 
density field! This is a dramatic illustration of the shaping 
of the cosmic web and the dominant role of voids in the dy- 
namics at these particular scales. The velocity field at scales 
of 4 h~^Mpc is highly coherent and characterized by lami- 
nar flows. Only at the boundaries of voids the velocity field 
becomes turbulent as we enter the dense regions dominated 
by non-linear interactions. 

Even more striking is the emergence of a hierarchy of 
structures in the velocity field as we probe smaller scales 
(right bottom panel of Figure [8]) . With the residual velocity 
field at scale of 1 h~^ Mpc the large central void seen as a 
coherent expanding domain at scale of 4 /i~^ Mpc fragments 
into several smaller voids. The nesting character of the ve- 
locity field is (as fas as the authors know) shown here in an 
explicit way for the first time. The expanding sub-domains 
embedded in the larger ones correspond to sub- voids in the 
density field. By using the residual velocity instead of the 
original velocities we are able to study the contribution of 
different scales in the cosmic web corresponding to different 
levels in the hierarchy of cosmic structures. Note that the 
parent void is perfectly delineated by its sub-voids except 
close to the boundaries where non-linear interactions make 
the velocity field turbulent. In the present work we provide a 
qualitative description of the laminar or turbulent character 
of the velocity field. This will be addressed quantitatively 
and in more detail in a future paper. 

The velocity field at the boundaries of the large voids 
tends to be highly turbulent refiecting the non-linear nature 
of the dense filaments and clusters surrounding the voids. 
This can be better seen in the residual velocity field at scale 
1 h~^Mpc (Figure [SI bottom right panel). The level of tur- 
bulence in the velocity field can be almost predicted by the 
local density (Figure (8] top left panel). This is not surprising 
since the high-density regions are dominated by small-range 
non-linear interactions. What is more interesting is the ve- 
locity field at the boundaries of the structures delineated by 
the sub and sub-sub- voids. The velocity field at the interior 
of the voids is practically laminar everywhere with the only 
exception of the boundaries delineating the voids at the top 
of the hierarchy. The highly coherent velocity field inside 
voids even at small scales refiects the fact that the internal 
substructure of voids is still close to the linear regime and 
there has been little shell crossing at the boundaries of the 
sub and sub-sub-voids. 

The magnitude of the residual velocity field inside voids 
at scales of ^-^ 1 h~^Mpc (Figure [9] right panel) is of the order 
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Figure 8. Hierarchical velocities inside voids. The top left panel shows a slice of the density field from which the central area is selected 
(yellow square). The right panels show the structures in the velocity field using the particle advection technique (see text for details) for 
three cases: original velocity field (right top) and residual velocity field at scales 4 and 2 /i~-'^Mpc (middle top and right bottom panels 
respectively). Left panels show again the velocity field for 4 and 2 /i~-'^Mpc and superimposed the Cosmic Spine showing the top level as 
a thick black contour (middle left panel) and the top+bottom levels as thick and thin black contours respectively (bottom left panel). 
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Figure 9. Magnitude of the velocity field for the original velocity field (left panel), and the residual velocity field (see text for details) 
at scales 4 (center panel) and 1 h~^Mpc (right panel). The magnitude corresponds to the full three-dimensional componets in contrast 
to Figure [8] where we only use the x,y components. The color bars on the top show the full dynamic range until the saturation point 
indicated by the highest value of the color bar. The lower color bar shows 1/4 of the full dynamical range. The matter distribution 
delineated by white filled isodensity contours at 6 = 0. The black square shows the region corresponding to Figures [2l [3l [SlfTO] and \TT\ 
Note tha the saturation point is 600, 200 and 100 km/s for the left, center and right panels respectively. 



of 10-20 km/s while at the interior of filaments at the same 
scales one can easily find velocities of > 100 km/s. This is 
a fundamental difference between voids and other environ- 
ments. Inside voids at small scales the local dynamics are 
still linear and strongly correlated to the underlying density 
field even at small scales. This may have important reper- 
cussions to the mass accretion in void galaxies as we will see 
in coming sections. 



6.2 Velocity field and the Cosmic Spine 

The expanding domains in Figure [8] have some interesting 
properties. They form a contiguous network of well defined, 
self-contained, locally expanding domains. There are no iso- 
lated structures in the velocity field just as in the density 
field i.e. one does not find a disjoint filament inside a void 
neither in the velocity or in the density field. Even if a 
structure appears not to be connected to anything, a higher- 
resolution inspection will always reveal it to be connected by 
a tenuous web. There is also no "leaking" of matter between 
adjacent voids. Once matter reaches the void boundary it 
stays there as it continues its evolution towards walls, fila- 
ments and finally clusters. The picture that emerges from 
the velocity patterns in Figure [8] is that of a network of con- 
tiguous and closed cells. The cellular nature of the Cosmic 
Web is perhaps the main reason behind the intimate rela- 
tion between the geometry and topology of the density field 
and why topological methods such as the Watershed Void 
Finder and the Cosmic Spine are able to identify geomet- 
rical structures in the density field by their topology. The 
collapse of structures described by Zel'dovich generates a 
range of geometries but it is their context within the Cosmic 
W eb that defines the ir connectivity and topology as hinted 
bv lBond et aP (|l996l l. In Sections [T2l and [STHI we discussed 
the motivation behind the use of a hierarchical approach to 
the identification of structures. There we argued that the 



over segmentation problem in the watershed void finder was 
not an intrinsic problem but (in our particular cosmologi- 
cal context) a direct consequence of the hierarchical nature 
of the cosmic web. It is then natural to expect that every 
feature in the density field, provided the density field was 
properly reconstructed, should be identified to a real cosmic 
structure. Figure [8] shows this very clearly. The top level of 
the Cosmic Spine hierarchy is shown on top of the residual 
velocity field and 4 h~^Mpc (middle right panel). There is 
an almost perfect match between the structures in the veloc- 
ity field (and density field, see Figure [2]). One must keep in 
mind that the Cosmic Spine method makes no assumption 
on the geometry of the density field and only relies on the 
topology of the local minima. Its success in the identification 
of structure relies entirely on the close relation between the 
geometry and topology of the cosmic web. If the cosmic web 
had different topological properties it would be impossible 
for the Cosmic Spine to identify any feature, for instance 
any disjoint structure (e.g. an isolated filament /wall) would 
not produce any change in the topology of the density field 
and would not be detected. 



A careful inspection of the middle right panel of Figure 
[8] shows that there are some voids that are not identified 
by the Cosmic Spine. These are structures that correspond 
to a different range in scales than the one covered by this 
particular hierarchical level. The bottom panel of Figure [8] 
shows the residual velocity field at scale 1 h~^Mpc and the 
top and bottom levels of the Cosmic Spine. Now the sub- 
domains are better correlated with the bottom level because 
we are in a different range of scales. Even at this scales the 
top level Cosmic Spine is closely delineated by the much 
smaller sub- voids. 
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Figure 10. Density and velocity profiles in the interior of a void. Top panels show a slice of the density field across the simulation box 
and the structures in the residual velocity field structures at scale 1 /i~-'^Mpc visualized with the particle advection technique (top left 
and right panels respectively). The top level of the hierarchical spine is superimposed on the residual velocity field. A vertical sampling 
line crosses the central void from point a to b. The sampled magnitude of the residual velocity field is shown in the middle right panel 
where the magnitude of the velocity field is decomposed into its residual components at scales 8,4, 2 and 1 /i~-'^Mpc (dotted dashed, solid 
black and solid red respectively). A small rectangle surrounding the vertical sampling line in the residual velocity field is highlighted in 
the top right panel and also shown in the middle right panel as the background for the residual velocity profiles in order to show the 
correspondence between the ID profile at 1 /i~-'^Mpc (solid red line) and the 2D structures in the velocity field at the same scales. The 
bottom left panel shows the density field along the vertical sampling line. 
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6.3 Residual velocity field profiles 

A more direct comparison between the density and velocity 
fields is presented in Figure [10] where we show the residual 
velocity profile at scales 8, 4, 2 and 1 h~^Mpc. As noted be- 
fore the structures in the velocity field seem to trace the 
hierarchy of voids more clearly than the density field by ex- 
plicitly showing their expanding character. Note the large 
variation between the different scales in the residual veloc- 
ity field. The large scales have larger amplitude as expected 
from the velocity power spectrum. 

The selected void has a density of ^ ^ —0.9 at its center 
and its density remains relatively constant until the regions 
close to its boundaries. The velocity field on the other hand 
shows a strong variation as we cross the void with a clear 
outflow pattern. This behaviour is seen only at scales > 
4 h~^Mpc (dotted, dashed and black solid lines in middle 
panel of Figure fTO)) At small scales the residual velocity field 
has a similar behavior than the density field. 

The residual velocity profile gives more information 
than a simple velocity profile where all the small-scale low- 
magnitude details would be lost in the high magnitude of 
the raw velocity field. Note the correspondence in the mag- 
nitude of the velocity profile at small scales (solid red line) 
and the 2D velocity structures in the top right panel. The 
solid red line crosses the y axis at the origin and bound- 
aries of the expanding domains in the background image. 
The same structures can be also seen in the density field 
(bottom right panel). 



7 VOID GALAXIES 

In the previous section we have seen that the velocity field 
inside voids is a superposition of hierarchically expanding 
domains dominated by large scales flows. Now we focus on 
the small scale velocity field that defines the local dynamics 
around the low-mass haloes found inside voids. 

Figure [11] shows a slice of the density field (top left 
panel) and the FoF haloes identified inside a slice of 2 
h~^Mpc of thickness. The haloes are arbitrarily scaled with 
their virial radius. Most of the haloes in the slice are located 
in filaments and clusters. There is a strong correlation be- 
tween the mass/size of the halo and its local density. The 
most massive haloes invariable occupy the dense regions and 
completely avoid the interior of voids. Instead voids are pop- 
ulated by low-mass haloes as one would expect given their 
underdense environment. This low-mass population is the 
combined result of the low-density environment which gives 
a slower gravitational development and the super-Hubble 
local dynamics around the haloes that decreases the radius 
from which the halo can accrete mass flAragon-Calvo et al.l 

In what follows we focus on the central void in Figure [TT] 
from which we select a halo far from its boundaries (the halo 
is shown in the top left panel enclosed in a yellow circle). 
The void halo has a mass of 2.79 x 10^° h~^ M© and a radius 
of 89 /i~^kpc. Figure [12] shows the density inside increasing 
spheres center in the void halo. At very small radius the en- 
closed density profile is dominated by the compact halo but 
at r ^ 1 h~^Mpc it rapidly falls to the mean void density 
at 6 ^ —0.9. After this point the profile slowly increases as 



the denser nearby walls and filaments are included in the 
sampling volume. This profile does not however shows the 
"real" profile of the void because the halo is not located at 
its center. At a distance from the halo of 8 h~^Mpc the en- 
closed density reaches again the mean density and continues 
increasing until it reaches its peak at K^ 9 h~^Mpc. If we 
ignore the contribution of the halo at small radius this would 
be a textbook void profile. This halo is a classical example of 
the cloud-in-void scenari o. Compare Figure [T2l with t he sec - 
ond panel in Figure 6 of ISheth fc van de WevgaertI (|2004l ). 
Both curves have a similar behavior even though our profile 
was computed from the present-time density field. 

7.1 Matter accretion in void haloes 

The dynamical environment around haloes is shown starting 
at the top right panel in Figure [TT] and continuing clockwise 
for the residual velocity field at scales 4,2 and 1 h~^Mpc 
respectively. At scales 4 and 2 h~^Mpc the velocity field 
around the void halo is basically laminar and the halo is 
simply dragged away from the void in the bulk flow. These 
flows have magnitudes of the order of 50-100 km/s (see Fig- 
ure [9|). 

At smaller scales the velocity field sudlendly presents 
a dramatic change. Instead of a large-scale laminar flow we 
now have small-scale coherent streams emanating from the 
centers of very small adjacent voids and reaching the halo in 
an anisotropic pattern. The small voids around the halo de- 
fine a web of tenuous filaments with the halo sitting at their 
intersection. This structure can be seen in both the den- 
sity and velocity field but the velocity field explicitly shows 
how the matter is locally shaped into this anisotropic con- 
figuration. This low mass halo is the dominant component 
in the local dynamics in a similar way as a cluster with a 
mass several orders of magnitude larger would be at much 
larger scales. The halo also accretes mass from its connected 
tenuous filaments again resembling the accretion of matter 
in massive clusters. Inside the void clean streams of matter 
feed the halo at very small rates, the residual velocity field 
at 1 h~^Mpc has a magnitude of the order of 10-20 km/s 
in contrast to the dense megapar sec- long filaments where it 
can be as high as 100 km/s at same scales (See Figure [9]). 

The most striking feature of the velocity field around 
void galaxies is its highly anisotropic and coherent nature. 
The local mini web around void haloes marks a fundamen- 
tal difference with similar haloes in denser environments. 
Here the halo dominates the local dynamics allowing it to 
feed from a well defined tenuous web. A similar halo located 
inside a filament will be embedded inside the larger dense fil- 
ament with non-linear turbulent dynamics. The high coher- 
ence and low magnitude of the velocity field around void ha- 
los may help explain the sustaining of the delicate extended 
gas disks characteristic of "isolated" galaxies (Kreckel et al.l 
l201lh and the also tenuous "rings" of gas observed around 
galaxies in low-density environments ( Stanonik et al.ll2009l V 

While the particle resolution of our simulation does not 
allows us to make accurate measurements of mass accretion 
we can make a first order estimation for this particular halo 
based on the local density and velocity field. The local web 
around the halo drags matter from its surroundings effec- 
tively increasing the local density. From Figure [12] we can 
safely assume a local density of 0.5 the mean density of the 
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Figure 11. Dark matter haloes superimposed on the density and velocity fields. The let top panel shows a slice of the density field and 
on top of it all the FoF haloes closer than 1 /i~-'^Mpc from the slice. We select a halo inside the central void with a yellow circle (top 
left panel) for subsequent analysis (see text for details) The residual velocity field is shown at scales 4,2 and 1 /i~-'^Mpc in the top right, 
bottom left and bottom right respectively. We show the same FoF haloes as in the density field slice. The velocity field was visualized 
using the particle advection technique (see text for details). 



Universe in the immediate neighborhood of the halo. The 
velocity field around the halo is of the order of ^ 20 km/s 
which gives a gas accretion rate inside a sphere of 1 /i~^Mpc 



provides a highly efficient me chanism to inject cold ^as into 
the i nner region of the halo (JKeres et al.1 l2005l : IPekel et alj 
l2QQ9l ). 



ion fraction of 0.04. If we assume that most of the accreted 
gas is converted into stars then we have similar values for 
star formation rates th an those obtained in high-resolution 
simulated void haloes ([Kreckel et al.ll201ll ). Perhaps more 
important than the magnitude of the accretion rate is the 
steady and coherent nature of the accretion process which 



8 CONCLUSIONS AND DISCUSSION 

In this paper we introduced the most recent implementa- 
tion of the Hierarchical Cosmic Spine method and applied 
it to a cosmological N-body simulation in order to decom- 
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Figure 12. Enclosed density profile centered in the halo selected 
in Figure [TT] 



pose the density field into its basic geometrical components. 
We focused our analysis on voids identified on a three-level 
hierarchical space. 

The size distribution of voids has a log- normal behay - 
ior in agreement with previous findings (jPlaten et al.ll2007h . 
Here we went one step further and computed the void size 
distribution for the three hierarchical levels computed in this 
work. At the top of the hierarchy the peak of the void size 
distribution is ^-^ 8 h~^Mpc which is very similar to estimates 
of th e radius of the local void (Nasonova & Karach entsevI 
l201lh indicating that the local void is average in size. The 
size distribution of sub- voids and sub-sub- voids have peaks 
at smaller scales as expected from their hierarchical nature. 
The void size distribution suggest that these structures are 
delineated by galaxies with luminosities in the ra nge -16 < 
Mr < -20 according to (JFurlanetto fc PiranllJOOeh . 

Voids are structures with very complex internal dynam- 
ics. By removing the large-scale component of the velocity 
field we illustrate how the hierarchy of voids in the matter 
distribution is reflected in a hierarchy of expanding domains 
in the velocity field. The notion of voids as simple expanding 
domains is far from accurate. Instead, voids have a rich in- 
ternal hierarchy of expanding sub-domains where the main 
expansion defining the parent void breaks down into smaller 
expanding sub-domains at smaller scales and so on. Further- 
more we also showed that the structures seen in the density 
and velocity field as well as their hierarchy are successfully 
identified by the Hierarchical Cosmic Spine method. Our 
method is unique in that it not only identifies structures in 
a hierarchical way but it also gives their explicit hierarchical 
relations. 

The velocity field inside voids can be characterized as a 
laminar flow dominated by large scales. The fluctuations in 
the residual velocity field at scales < 1 h~^ Mpc are one or- 
der of magnitude lower than the fluctuations at scales com- 
parable to the void size. The overall internal dynamics of 
voids are largely dominated by scales > 4 h~^Mpc while the 
smaller scales play a minor role. 

We presented an exploratory study of haloes inside 
voids. Contrary to voids haloes are not affected by the large 
scale motions. They simply move together with their envi- 



ronment in large-scale and large-amplitude bulk flows. Small 
scales on the other hand play a major role by defining the 
dynamics around the halo and the properties of the galaxies 
they host. We found that haloes inside voids feed through 
highly coherent laminar streams of matter in sharp contrast 
with the environment inside filaments and clusters where the 
velocity field is turbulent. These streams produce a highly 
anisotropic matter accretion pattern around the halo and 
may provide an efficient gas feeding mechanism with low 
accretion rates that can be sustained for large periods of 
time. This steady and coherent gas accretion may offer a 
natural explanation for some of the observed properties of 
galaxies in voids and low-density environments such as ex- 
tended gas disks, polar disks, higher gas content and star 
formation rates. I a forthcoming paper we will explore in 
more detail the effect of the small-scale velocity field in the 
gas accretion of galaxies in voids and how this affects their 
observed properties. 

Voids represent a unique opportunity to probe the 
structure and dynamics of the early Universe at small scales 
by doing local observations. Future ultra-deep surveys cov- 
ering the volume of nearby voids may one day unveil their 
rich substructur e and dynamics. In this direction the Void 
Galaxy Survey (JKreckel et al.l [2OIII I is leading the way in 
the study of this fascinating objects with already some tan- 
talizing results like the possible finding of the first inter- void 
mini filament in the Ursa Minor I void seen as a long tenu- 
ous bridge of matter con n ecting three galaxies (VGS-38 in 
Figure 7 of iKreckel etHI (|2Qllh ). 
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